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Abstract. We compute the semi-global symplectic invariants near the hy- 
perbolic equilibrium points of the Euler top. The Birkhoff normal form at the 
hyperbolic point is computed using Lie series. The actions near the hyper- 
bolic point are found using Frobenius expansion of its Picard-Fuchs equation. 
We show that the Birkhoff normal form can also be found by inverting the 
regular solution of the Picard-Fuchs equation. Composition of the singular 
action integral with the Birkhoff normal form gives the semi-global symplectic 
invariant. Finally we discuss the convergence of these invariants using both 
analytical and numerical arguments. 
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The Euler top is a classical Hamiltonian system which describes the rotational 
motion of a rigid body in free space in the absence of a potential. In the centre of 
mass frame this can be considered as the rotation of the rigid body about a fixed 
point. This system has been treated extensively in the literature for the last 300 
years, so we only give a few selected references here [lOi The Hamiltonian H 
can be written in terms of the moment of inertia tensor M := diag(0i, 02, Q3) and 
the angular momentum L := {Li, L2, L3) as 



Without loss of generality, we assume the ordering < Oi < 82 < 63. The 
principal moments of inertia satisfy triangle inequalities Oi < 82 + 63 and so forth 
cyclically. We exclude the (degenerate) cases 8i = O2, 82 = 83 since then there is 
no hyperbolic equilibrium. This way of writing the Hamiltonian is not symplectic as 
the angular momenta are not canonical variables. Instead the system has a Poisson 
structure with Poisson structure matrix 



MSG: 37J35, 37J15, 37J40, 70H06, 70H08, 37G20 



1. Introduction 



(1.1) 




/ -L3 



L2 





so that the differential equation is written as 



L = FbVH = i X VH 



Due to conservation of angular momentum, the dynamics of the Euler top takes 
place on spheres of constant magnitude of angular momentum ||Z|| — i. In fact 
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the total angular momentum is a Casimir C{L) := \\L\\'^ of the Poisson structure, 
PeVC = 0. The level set of the Cashnir 

Ce := {LeR^:C{L)^e^} 

is a sphere and the level set of the Hamiltonian, the energy surface 



|i e M^: H(i) = h| 



is an ellipsoid. The solution curves are given by the intersection CeHSy^. The (non- 
degenerate) intersections of the sphere with the ellipsoid give concentric ellipse-like 
curves, centred about elliptic equilibria of which there are four. There are two 
separatricies, the intersections of which occur at the two hyperbolic equilibria. In 
total there are six equilibria given by ±£Li, i £ {1,2,3} where Li is defined as 
the unit vector along the Li axis. We are interested in the unstable hyperbolic 
equilibria (i = 2), corresponding to a steady rotation about the principal axis of 
inertia corresponding to the middle moment of inertia 02. The eigenvalues of the 
linearisation about the hyperbolic equilibria are ±A, where 

(1.2) X = £^iA-B){B-C) , 

using the abbreviations 

(1.3a) A:^e^^ 
(1.3b) B := 

(1.3c) C :^ 63^^ 

for the reciprocals of the principal moments of inertia (hence O<8i<02<03=S> 
< C < B < yl). For later use we also define the dimensionless parameters 



B-C 
A-B 
1 



(1.4a) p := 

(1.4b) K := 

P 

so that p ^ \/p and k — — k when 81 and 83 are exchanged. The involution 
p -l/p leaves k invariant. Note that n"^ = p"^ + p-^ - 2 = {A~2B ^ CY/{A - 
B)/{B — C) is rational in the moments of inertia. When restricting to the physical 
range p > the equation k = p— \/ p yields the unique solution p = \ (k+Vh^ + 4 ). 

In this work the main aim is to calculate the semi-global symplectic invariants 
of the Euler top near the hyperbolic equilibrium point. In [S1|T5], Dufour, Molino 
and Toulet introduce the classification of integrable systems using their semi-global 
symplectic invariants. Their approach considers the triple (A^,a;, J-"), where is a 
two-dimensional manifold, w is a symplectic 2-form, and is a Morse foliation given 
by the levels of a Morse function F. Then the equivalence between two integrable 
Hamiltonian dynamical systems with one degree of freedom is introduced in [6J by 
the following definition: 

Definition 1. Two triples (A^ 1 , Wi , J^i ) and (AI27 ■^2) -^2) £^re said to be equivalent if 
there exists a symplectomorphism between A4i and AI2 that preserves the foliation. 

In formulas this means there is a symplectic diffeomorphism 4* : A^i — ^ A^2 
and a diffeomorphism S : K — >■ K such that Fi o 'if = S o F2 . They introduce the 
semi-global symplectic invariant as the regular part of the action integral near a 
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separatrix written in a particular canonical coordinates system defined using the 
Birkhoff normal form at the hyperbolic equilibrium point. The main theorem in [B] 
then is 

Theorem 2. In a full neighbourhood of a separatrix of the same topological type 
two systems are equivalent if and only if their semi-global symplectic invariants 
coincide. 

This theorem motivates us to calculate the semi-global symplectic invariant of 
the Euler top. This is the first step that may lead to a non-trivial equivalence 
between the Euler top and other Hamiltonian systems with one degree of freedom. 
We can, e.g., decide whether a particular Euler top (it is a three parameter family) 
is equivalent to say the pendulum, for which the invariants have been calculated in 

The theory of semiglobal symplectic invariants was further developed by San Vu 
Ngoc in [T^ who extended it to focus-focus equilibria for system of two degrees of 
freedom. The first explicit computation of semi-global invariants near a focus-focus 
point was done by DuUin in for the spherical pendulum. Hyperbolic-Hyperbolic 
points were treated in [7, where the C. Neumann system was used to illustrate the 
general theory. Other notions of equivalence have been used to study the Euler top 
in [21 [mill]. 

Plan of the paper: we will be calculating the Birkhoff normal form algorithmically 
using Lie series. For one degree of freedom, the Birkhoff normal form about a 
hyperbolic equilibrium is a series in powers of the regular action only. Next we 
show that the action integrals satisfy a Picard-Fuchs ODE. Then we solve this 
ODE using the method of Frobenius, and arrive at series expansions for the action 
integrals. We find that if we invert the regular action integral then we recover the 
Birkhoff normal form. Once we have the Birkhoff normal form, we calculate the 
symplectic invariant by composing the singular action integral with the Birkhoff 
normal form. 



2. Calculating the Birkhoff normal form 

In order to calculate the Birkhoff normal form of the Hamiltonian of the Euler 
top at the unstable equilibrium, we introduce local canonical variables {q,p) with 

1 



symplectic structure S" .— , ^ ^ 

Lemma 3. The Poisson map $ : M'^ — > defined by 

(2.1a) g = Arg(L2+iii) 

(2.1b) p = L3 

maps the Hamiltonian H{L) of ( |1.1| with value h and Poisson structure Pe on the 
symplectic leaf Ce into the standard symplectic structure S in canonical variables 
{q,p) and Hamiltonian H{q,p) with value h given by 

(2.2) H{q,p) = I [p' (C - fiq)) - £2 (B _ /(5))] 

where f{q) -.— A sin^ q + B cos^ q. 
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Proof. Substituting equations ( |2.l| ) into (2.2 I and using i'^ = L\ + L\ + L\ yields 



the required Haniiltonian up to the constant ^P'B which is absorbed in the 

definition of H for convenience, so that H{^{L)) + \i'^B = H(-E/). The hyperbolic 
equilibrium ILi is mapped to the origin ((7,p) — (0,0), while —IL2 is mapped to 
{q,p) = (tt, 0). To derive the new symplectic structure, compute the 2x3 Jacobian 
matrix and verify that D^Pe {D't)^ — S. □ 

Note that this transformation to "cylindrical coordinates for the sphere" is not 
defined globally on sphere, but only on the punctured sphere with the two points 
±££3 (where Li — L2 = 0) removed. However, the transformation is valid near 
the unstable equilibria ±£L2 and in a full neighbourhood of their separatrix. An 
alternative transformation introduces Li as momentum instead of L3, and the cor- 
responding Haniiltonian H{q,p) has A and C interchanged. 

Lemma 4. Using 1/A as units of time, I as units of angular momentum, andl/X 
as units of moment of inertia, the Hamiltonian in non-dimensional form is 

(2.3) H{q,p) = i {-p\p + sin^ q) + p-^ sin^ q) . 

The proof is a simple calculation. The new scaled angular momentum p = pjl 
is dimensionless, as is the value of the Hamiltonian h = h/{X£). 

Remark 5. Each time a transformation is done, the variables change. For clarity 
and simplicity of notation we use the same letters {q,p) for old and new variables, 
but it should be noted that each transformation introduces different variables. In 
our notation the tilde designates quantities with dimensions, while from this point 
onwards we use non-dimensionalized quantities H, q, h without tilde. Sans-serif font 
is used to designate quantities in the original Poisson system, so that h — -^£^3 = 
h = hX£. 

Remark 6. The area form on the original sphere in these variables is £dq A dp — 
i'^ dq A dp so that the scaled symplectic area dq A dp differs from the true area on 
the Casimir sphere Ci by a factor of £^ . 

The original Hamiltonian has a group of discrete symmetries generated by L,; — >■ 
—Li, i — 1,2, 3. In the canonical variables these correspond to g — >■ —q, q ^ tt — q, 
p -> —p, respectively. The global analysis we are going to present later is simplest 
if there is only a single hyperbolic equilibrium on the separatrix in question, and 
therefore we are going to consider the Euler top modulo its discrete symmetry group. 
Any two pairs of the three discreet symmetries generate the group of symplectic 
discrete symmetries of the Euler top, which is isomorphic to Klein's Vierergruppe 
V = Z2<S)'Z2- In the canonical variables a possible choice of generators is Si{q,p) = 
{~q, —p) and S2{q,p) — (7r-|-g,p) which are both involutions. A fundamental region 
for the quotient of the cylindrical {q,p) phase space [— tt, tt) x (—1, 1) by the group V 
generated by and S2 can be chosen as the positive quadrant {q,p) G [0, tt) x [0, 1). 
This corresponds to a quarter of the original sphere 

V (Ce) := {L eCi-. < Li < £,-£ < L2 < £,0 < L3 < £} . 

We Taylor expand the non-dimensional Hamiltonian H{q,p) about the origin 
{q,p) — (0,0) for analysis near the equilibrium L — £L2. The quadratic terms are 
\{—p^p + q^p^^)- The Williamson (linear) normal form of the hyperbolic equilibria 
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is found after a symplectic linear transformation (e.g. as outlined in [3 ). Although 
the Williamson normal form is unique up to the overall sign of the qp term which 
we chose to be positive, the transformation is not; we chose to perform a symplectic 
scaling q — ^J~pq^]3 vl \fp followed by a rotation by — 7r/4, so that the positive 
quadrant in the new coordinates corresponds to positive Hamiltonian. With this 
convention the Williamson normal form becomes unique. 

Lemma 7. The symplectic linear transformation 

gives the Williamson's normal form H^, of the Hamiltonian H as 

H.=qp- - q^f -^^{p + qY + 0(6) 

In order to arrive at the BirkhofF normal form, we use the method of Lie trans- 
forms to remove terms that are not powers of qp. We now state the Lemma that 
gives the recursive algorithm as given in |I12| : 

Lemma 8. The Lie transform generated by W transform,s Williamson's normal 
form into the Birkhoff normal form H* using the double-indexed array {^j}^ -^^^ 
related recursively via 



where {•,•} is the standard Poisson bracket in the canonical co-ordinates {q,p). 

The generating functions Wi are chosen to be polynomials in g,p of degree 2(^+1) 
with coefficients chosen in such a way that they remove terms that are not powers 
of qp. The algorithm is implemented in Mathematica; we state the main results in 
the following Theorem: 

Theorem 9. The Birkhoff normal form up to order 7 of the Euler top at the 
hyperbolic equilibrium point is given by 

^ '^--^'^ -^^'^ 12^^ 

3 {n^ + 4) (11^^ + 12) In (^^ + 4) (9^^ + 20) 
^ ' 1024 2048 

(^2 + 4) (527^4 + 1776^2 + 720) ^ , 

16384 ^ ^ 

where J = qp in the new variables. 

Note that the parameter dependence on the right hand side is only through 
the dimensionless parameter n = p — 1/ p and the power series is in terms of the 
dimensionless action J. 

3. Action integrals via Picard-Fuchs equation 

The essential step in the calculation of the semi-global symplectic invariants is the 
computation of the action integrals, which are given by complete elliptic integrals 
in the case of the Euler top. Since we are interested in the series expansions of these 
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(a) Casimir sphere Ci 



(b) Canonical phase space 



Figure 3.1. Areas enclosed by the same closed orbits after dis- 
crete symmetry reduction (A) on the quarter Casimir sphere 
L E V (Cf ) and (B) on the scaled and discrete symmetry reduced 
canonical phase space {q,p) G [0,7r) x [0,1). The parameter is 
K = 1/2. The thick black line is the separatrix with h = 0. The 
lighter (blue) shade below the separatrix shows the action of an 
orbit with h > and the darker (red) shade above the separatrix 
shows the action of an orbit with h < 0. The symplectic invariant 
is calculated from these areas in the singular limit /i — > 0. 

integrals the most natural approach is not through the integral itself, but instead 
through the so called Picard-Fuchs ODE that the integral satisfies. The derivation 
of the Picard-Fuchs equation proceeds in a way similar to [9^ . Frobenius expansions 
of this linear ODE then gives the desired series. This gives a basis for the vector 
space of solutions of the linear ODE, and in a second step the particular solutions 
corresponding to the action integrals of the Euler top are found. 

Lemma 10. The scaled and discrete symmetry reduced action of the Euler top with 
scaled energy h normalised to at the unstable equilibrium is a complete elliptic 
integral on the curve 



r := {{z,u) e C^\u^ 



{2h-z)w{zf, w{zf 



z{z^ + KZ~ 1)} 



given by 




along the cycles /3± for zth > as specified in Fig. \3.S\ 



Proof. We derive the action using the Hamiltonian H[q,p) in canonical variables 
{q,p). Solve H{q,p) — h for p and integrating with respect to q in order to find the 
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lighter (blue) area in Fig. 3.1 gives 



1 / 2h-p-^fim^q 



91 



where qi = 11/2 — q^^ and q2 = 7r/2 + go ^re the roots of the numerator of the 
integrand with go = cosr^{^2hp'). Upon performing the change of integration 
variable z = p^^sin [q) we arrive at the required integral over the differential 
1-form :— ^^^~)^ dz. The new denominator [w(z)]'^ ~ z (z^ + kz — l) = 

z(z + p){z — ^) depends only upon the single parameter k = p — Note that 
when q traverses the interval [91,92], then z completes the closed path /3+ once. 
The complex closed loop integral around /3+ gives twice the real integral along the 
interval [2h, l/p], and so altogether Ip^ = ^ (. 

A slightly more complicated argument applies in the case h < 0. To get the 



lighter shaded (red) area in Figure 3.1 the integrand is (1 — p) dq instead of just pdq 
and the real integration interval is q G [0,7r]. Complexifying this integral gives the 
loop so that — \ — ^ §p C- On the Riemann sphere we have /3* +/3++7oo = 
where is a loop around the pole at infinity. Due to the non-zero residue of C, at 
infinity the constant term cancels and we get = C- I— ' 

Remark 11. The unsealed action / as a function of the unsealed and unshifted 
energy h can be rewritten in the symmetric form 

1 r / 2h - /2? ' 

The scaled action / depends on h and t only through h ^ (h — ^B£'^)/ (X£), and 
on A, B, C only through k. In the transformed variable the roots z — A,B,C 
correspond to the roots z — — p, 0, ^ of w, respectively. 

3.1. Derivation of the Picard-Fuchs equation. We now derive the Picard- 
Fuchs ODE of I{h). C, is an Abelian differential living on the complex manifold F. 
By deRham cohomology theory, there must exist a relationship between derivatives 
of C on F. In fact, we follow the route of Clemens in [5J and find that there exists a 
linear combination of the derivatives of ^ that equals a total differential. However, 
we cannot exploit the same simplifications as Clemens does since our independent 
variable is fixed to be the energy h, because otherwise we would lose the connection 
to the Birkhoff normal form (see below). A similar approach was taken in |9j. The 
relation between the differentials is given in the following Lemma: 

Lemma 12. There exists a function v{z) meromorphic on F and coefficients Ci 
such that 

(3.1) = 



i=0 



Proof. Observe that 



5 ^ 

w(z)(2ft,- z)2 "S^Cj— ^ = \-coz^ + (6coh + cAz'^ + (-UcqH'^ - Acih + C2)z 



i=0 



-(8co/i^ + 4cift.^ - 2C2/1 + 3C3)] dz 
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(b) h<0 



Figure 3.2. The a and /3 cycles, choices of branch cuts and 
branch points of T in the z E C plane for /i > and h < 0. 



Now choose v{z) := — = ^, which is meroniorphic on T. and has 

differential 

w{z){2h - z) 2dw = [(3/i + iK)z^ + [2nh - \)z - h\ dz. 

Equating the coefficients of the polynomials in z and solving for yields the unique 
solution 

Co = 

ci — 3ft. + |k 

C2 = 12/(,2 + 4k/i - 1 

C3 = ft(4ft2 + 2/iK - 1) = ^ \w{2h)f 

and thus by construction we have proven the Lemma. □ 

Now we are ready to derive the linear and homogeneous Picard-Fuchs equation 
for the action 1(K): 
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Theorem 13. The scaled action I{h) satisfies the Picard-Fuchs equation 
(3.2) [w{2h)f I"'{h) + 2(12/1^ + 4k/i - 1)1" {h) + {Qh + k)I' {h) = 

with the scaled energy h as independent variable. 

Proof. To obtain the Picard-Fuchs ODE, perform a closed complex contour integral 
to both sides of equation (3.1). By definition we have = 47r/, and given v is 
meromorphic on T the residues of dv are vanishing, so that the right hand side gives 
Sdv — for any closed integration path. □ 



3.2. Solving the Picard-Fuchs equation. Clearly (3.2 1 is an ODE in /' thus 
/ = fca is a constant solution. To lower the order we introduce the scaled period 
T{h) = 2TrI'{h), which has the first kind differential dz/u on F. The scaled period 
T{h) hence satisfies the second order linear homogeneous ODE 

(3.3) T"(h) + 2 5 r(h) + —^T{h) = 0. 

[w{2h)f [w{2h)f 

It is interesting to observe that the leading coefficient C3 is proportional to [w{2h)f' . 

and thus in normalising the ODE the roots of w{2h) given by 2h e |o, — p, ^| 

become the (regular) singular points of the Picard-Fuchs ODE. Thus the partial 
fraction decomposition of the coefficient of T' simply is 




We are interested in series solutions at the singular point corresponding to the 
unstable equilibrium, namely h = Q. The general theory and procedure for solving 
(3.3 1 via the method of Frobenius can, e.g., be found in We seek series solutions 
of the form 

CO 

^a„(g)/i»+e 

ri=0 

where g is a root of the indicial equation. At the finite singular points the indicial 
equation is gi^ = 0. 

Remark 14. We find that the Picard-Fuchs ODE has a regular singular point at 
h = 00, confirming it is of Fuchsian type. The indicial equation for the singular 
point at infinity is 

qI^-2q^ + 1=Q 
with roots Qrx, = 51 | differing by an integer. 

By substitution of the Frobenius series into the ODE, the recursion relation for 
a„(p) is found to be 

(3.4) a„(g) = ^"^^ ^ (^(2^^ + - ^)an-i{Q) + (2n + 2q- 3)a„_2(f?)) 
Setting Q — Q and a„ :— a„(0) yields the Frobenius expansion of the regular solution 



00 



Tr{h) ^^anh" 

n=0 
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with coefficients obtained from (3.4) at g = as 

(3.5) a„ = (2 '^^'^ ~ + (2" - 3)a„_2j • 

We solve this second order recursion relation for a„. Without loss of generality, 
normalise the initial condition gq 1, and we require that a_i := 0. Thus we find 
that the next few coefficients are 



ai 



K 

2' 



«2 = ^(3«^+4), 
5 

as = — a«(5k2 + 12), 
35 

a4 = —(35^^4 + 120^2+48), 
63k 

as = — — (63^-* + 280^2 + 240) 
2048 

Theorem 15. The recursion for an is solved by 

1 f2n\ / 2n-2fc \ /K\"-2fc 



4" V n / -^-^ \k,n — k.n — 2k 



(i 



/c=0 

where ( . "j.) = iijffcT n = i + j + k is the trinomial coefficient. 

The proof of this theorem will be given later as a special case of the (approximate) 
solution of the more general recursion for a„(gi). 

Remark 16. It is interesting to note that the sum a„ can be summed to the hyper- 
geometric function 2F1 

{2n)lY „ /n-l n 2n - 1 4 



(n!)2 ; ^ ' \ 2 ' 2' 2 ' 

which is always terminating because n is an integer. It is well known that the 
complete elliptic integral of first kind T can be expressed in terms of the hypergeo- 
metric function as the function of the modulus of the elliptic curve, while here the 
coefficients of the Taylor series of T,.(/i) are given by the hypergeometric function. 

Since we have repeated indicial roots, we expect the second independent solution 
to be singular. The general theory (see for example [1]) says that the singular 
solution is of the form 

00 

T,{h) :=r,(/i)log/i + ^6„/i" 

n=l 

The recursion relation for the b„ at p = is thus given by 



where &„ := '^"^'^^^•^ 



(3.6) 



g=Q 



_{2n- 1) (2Ka„_i + Kn{2n - l)&„_i + 2n(2n - 3)fe„_2) + 4(4n - 3)a„_2 
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Along with the initial conditions on the a„, we also impose that b-i := and 
&o 0. The first few coefficients are 

bi = K, 

b, = l(21«2 + 20), 

^3 = ?^'«(l85«:2 + 372), 
96 

64 = TT^T (iSeSSK'^ + 56760^2 + 18672) , 

6144 ^ ' 

65 = (102501k'* + 416360^2 + 313680) . 

20480 

We were not able to find an explicit solution for an{g)- However, since we only 
need the derivative for an{g) at g = it is enough to find an approximate solution 
dn{g) that is valid up to terms of 0{g'^). Using this we obtain an explicit formula 
for 6„ = a' JO) = aj^)- 



Lemma 17. The recursion for a„(gi) given by (3.4| is approximately solved by 

~ f ^ 0« „ ,(£'+l/2)„ ig+l/'2.)n-k -2k 

(f?+l)2 {n-2k)\k\ 
to leading order in g, where (a;)„ is the Pochhammer symbol. 

Note that this formula reduces to the explicit formula for a„ given earlier when 
£1 = using the identity (l/2)„ = (2n - l)!/(n - l)!/22"-i. Hence the following 
proof will also prove Theorem [T5| 

Proof. We will show that a„(p) — ajg) + 0{g^) by induction. From the explicit 
recursion we find 

(l + 2g)^ 

and 

(l + 2g)(3 + 2g) (l + 2g)2(3 + 2g)2 2 
(2 + g)2 + 4(l + g)2(2 + g)2 ■ 
One can check that ai(g) = ai(g) and that the coefficients of in 02 (g) and 02 (g) 
coincide and that the constant coefficients of 02 (g) and 02 (p) are both 

| + |g + 0(g2). 

Now assume that the identity a„(g) = a,i(g) + ©(g^) holds for all 1 < n < m — 1 
for some fixed 2 < m e N. Thus we aim to show that under these assumptions, 
o-mig) — flm((?) + ^{g^) holds. We substitute the formulas am_i(g) and am-2{g) 
into the right hand side of (3.4 1 with n = m and need to verify that am{g) so 
obtained is equal to am{g) + ©(g^). 

Define := am(g)(2m + 2g + 1) so that 

, J ,0 J X 1 (2m + 2g+l)(2m + 2g-l) 

dm = gm(Kdm~i 2 + dm-2), whcrc g,n := ; — . 

[m + gy 

From a we obtain d which cam be written as 

L™/2J 

dm = Gm{g) Sm,m-2k{g)l^™ 
fc=0 
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where 

2"m!(g+ 1/2)^+1 _ 2(g+l/2)(,+,)/2 

The claim is dm = 9m{i^dm-i/'^ + <^m-2) + Inserting the series for d and 

collecting powers of k gives 

,m—2k ,m— 2fe— 1 

Dividing by Gm-i{Q)Sm-2,m-2k{Q) and using the identities 

Gm{g) _ 2m(m + g + 1/2) _ m 
Gm-iie) ~ (^ + rn)2 ~ 2rn + 2g- 1^" 

and 

Sm,m-2k{Q) _ m - fc + g - 1/2 Sm-l,m-2k-l{B) _ m - 2k 
Sm-2,m-2k{6) k Sm-2,m-2k is) k 

gives after some simplification 

g2 = im-l){m + gr 

{2m + 2Q+l)Gm-l{Q)Sm-2 ,m—2k 

Evaluating the factor of ©(g^) at g = gives 

(m—X)m} 7r(TO — l)TOm!(fc — l)!(m — 2fc)! 

~ (2m + l)G„_i(0)5„_2,rr.-2fe(0) " " 2™+ir (m + I) r (m - fc - i) ' 
which is non-zero for all values of m > 2. 

□ 

With the previous Lemma it is now straightforward to find an explicit formula 

for the coefficients &„ . 



Theorem 18. The recursion relation for bn is solved by 
bn = 



1 /2ny^Y 2n-2fc \ / ny-^k 

fe=0 



4" V n 
where 

fn,k ■= 20n + 20„_fe - 2H„ 

where H„ is the Harmonic number and 0„ its odd cousin defined by 

n ^1 

H„:=^l/fc, 0„:=^-— . 
fe=i fe=i '^'^ ^ 

Proof. Using the previous Lemma wc can simply differentiate an{g) and evaluate at 
in order to get 6„. The derivative of the Pochhammer function is given in terms 
of the digamma function ip, which for integer and half- integer values n can be 
expressed in terms of the Euler-Mascheroni constant 7 and the Harmonic numbers, 
tp{n) = —7 + H„_i, where H„ = J2k=i ■ Using the recursion H„ = H„_i + 1/n 
and H(l/2) = 2 — 2 log 2 the Harmonic number is thus also defined for half integers; 
explicitly = 2D„ +H_i/2- Denote by a^(g) the coefficient of k"~'^'' in a„(g). 

The logarithmic derivative of a^(g) at g = can thus be found as 

-frd = = ""-1/2 + Hn-fe-1/2 - 2H„ + 41og2 = 20„ + 20„_fe - 2H„ 

and this determines the "correction factor" for the coefficient of fi;^"-^ in 6„. □ 
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To obtain the solutions Ir{h) and Is{h) of the Picard-Fuchs equation we integrate 
Tj.{h) and Tg(h) term- by-term, respectively, and get 

(3.7a) Irih) := -L f T^^h) dh ^ ^h-+' 

2n I ZTT ^-^ n+1 

■' n=0 

(3.7b) h(h) ^ jTs{h)dh 

oo . 

(3.7c) = Irih)lo,h+-Y^ — [l 



n + 1 



where the integration constants are fixed by the requirements that /^(O) — and 
Is{h) ^ as /i ^ 0. 



3.3. Particular action integrals. Since (3.2 1 is a linear third order equation, 
there must be three linearly independent solutions. They are: the regular, singular, 
and constant solutions. Thus the general solution is an arbitrary linear combination 
of these, namely 

(3.8) I{h)=kilr{h)+k2ls{h)+k3, 

and upon differentiating 

T{h) = kiTr{h) + k2Ts{h). 
We seek to find the ki that give the particular solutions corresponding to the 



closed loops integrals along the paths (3± as specified in Figure 3.2 The ex- 
pansions obtained are normalised such that Tr = 1 + 0{h), 2'Klr ^ h + 0{h'^), 
Ts — log ft, + 0{h), 2ttIs = hlogh + 0{h), and so the leading terms are 27Tl(h) = 
2Trk3 + k2h\ogh + kih + 0{h'^) and T{h) = ki + k2\ogh + 0{h). Thus the constant 
fca is given by /(O), while both ki and k2 are given by the leading order behaviour 
of T{h). When T{h) is finite for ft —> then k2 — and fci is determined by r(0). 
Otherwise the leading order logarithmically diverging term and the constant term 
of T(ft) for small h determine ki and k2. 

The four particular solutions to be found are given by the integrals Ip^ and Tp^ . 
In order to find the correct linear combinations we need to evaluate these integrals 
in the limit ft — )■ 0. 

The f3 integrals at ft = are computed as real integrals. The f3 cycles are shrunk 
down, so that we are integrating along the real intervals 11 [Z^+J := (2ft,-) and 



p 

n[/3_] := (— p, 2ft). /^±(0) is an elementary and finite real integral that gives 

1 / dz 

47r fp^ y/l — KZ — z^' 



If dx 1 1 / ici\ 

^ / , = - tan^i (p=Fi) 

27r Jn\8+] J(x + n)(l/p-x) tt 



2^ Jn[p±] y/ix + p){l/p-x)' 
For the singular integral T^^(ft), the asymptotic behaviour for small ft is 



2w(z)V2ft - z' 7n[/3±] ^/x(x^^^2hj^ \/r 
1 
2 



1 / 64 
±log(±ft)T;^log^— ) +0(ft) 
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This can be found as follows. Define (^(x) := [I — kx — x'^] 2 = [[x + p){l/ p — x)] 2 
and split the integrals up as 

T,JH)=[ ^^P^ax^f ^ ,d.. 

Jn[p±] ^Jx{x~2h) Jii[p±] y/ x{x — 2h) 

The first integral is a convergent elliptic integral, and when ft, = it becomes 
elementary. For /3+ it gives — log ^ ^^^^f ^ , and for /?_ it gives log ^^^^5^^ The 
second integral is divergent when /i — >■ but elementary and can be integrated using 
hyperbolic trigonometric substitutions. For /3+ it gives — log + \ogh + 0{h), 

and for /3_ it gives log (2p) — log(— ft) + 0{h). Adding the two integrals gives the 
stated result. 

From these 4 integrals the coefficients ki can be determined as described above 
and the result is 

(3.9) = tI log (;^) Ir±Is + l tan-i {p^') . 

Remark 19. Notice that the actions related to approaching the separatrix from 
either side satisfy Ijs^ (0) + Ip_ (0) = ^ , corresponding to the fact that the total area 
of the symmetry reduced scaled phase space is tt. This area can also be seen as the 
residue at infinity of /(ft) 

27r (/^^ (0) + Ip_ (0)) = i • 27rilte {((ft)} = n. 

As mentioned earlier the actual area (without discrete symmetry reduction) en- 
closed by a single connected contour of H is twice as large, and each connected 
component appear twice. Undoing the scaling then gives 47r£^ which is the area of 
the sphere Ce which is the (5*0(3) reduced) phase space of the Euler top. 

4. The symplectic invariants 

Equipped with the Frobenius series expansions of the action integrals obtained 
from the Picard-Fuchs ODE, we can now calculate the semi-global symplectic in- 
variants of the Euler top. 

4.1. Revisiting the BirkhofF normal form. The Birkhoff normal form is a series 



for ft (J). From (3.7a I we instead have a series for the regular action in terms of the 
energy 

27r/^(ft) ^ h + -h^ + — (3k^ + 4:) + . . . 
4 16 

The regular action can be obtained by integrating ( over the a cycles 
(4.1) Uh) ^ £ CW = f'+J^'^ + Y^ (3'*' + 4) ft-V . . . 

where the series can be obtained by Taylor expansion and taking residues. We can 
omit the subscript ± for a since the two cases yield the same series expansion. The 
relation between the action integral and the regular Frobenius series thus is 

/„(ft) = 27r/^(ft). 

Thus by inverting this series we recover the Birkhoff normal form, so that we can 
identify — J. The fact that the Birkhoff normal form at a hyperbolic point is 
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given by the integral of the a. cycles is a general phenomenon, see [8] for a general 
proof. 

4.2. The semi-global symplectic invariant. The semi-global symplectic invari- 
ant cr(J) is the power series given by the regular part of the composition of the 
singular action integral with the inverse of the regular action integral. More pre- 
cisely there are two cases, depending on the sign of h: 

(4.2) 27r (/^^ o /-I) ( J) = ^± ± J log (± J) T J T a( J) 

Here A± is the area enclosed by the separatrix after discrete symmetry reduction, 
A± = 27r/^^(0) = 2tan"i {p^^) so that + A- = tt. Thus we obtain the 
semi-global symplectic invariant a: 

Theorem 20. The semi-global symplectic invariant up to order 7 of the Euler top 
with distinct moments of inertia reduced by discrete symmetry near the hyperbolic 
equilibrium is given by 

1, / 64 \ ^ 3k. 15k2+32. 5k (Uk,^ + 36) . 

^ ' 2 ^\k? + A) 8 96 512 

945k'^ -I- 4200^2 -I- 2672 . 7k (527k4 + 2960k2 + 3600) g 

10240 40960 
65709kS + 44604Gk4 + 80136Gk2 + 241664 ^7 , 

r + o (r) 

688128 ^ 



Proof. The action Ip{h) is given as the series expansions (3.9 1, (3.7 1 whose coeffi. 



cients where obtained from the Frobenius expansion of the Picard-Fuchs equation. 



The action Ia{h) is similarly given by (4.1), and the inverse of this series is the 



BirkhofF normal form H*{J). Composing Ip{h) with H*{J) gives the series of the 
action Ip in terms of the normal form action J, from which the symplectic invariant 



<^e{J) can be read off using the definition in (4.2 1. □ 



Notice that for k = (i.e. p = 1, i.e. A — B = B — C) the invariant is an 
even function of h, so that both sides of the separatrix, for positive and negative 
h OT J are the same. The linear term has maximal value of log 4 at k = 0. For 
positive K all higher order terms are negative. The series expansion of as has 
been numerically verified and agrees well with the values obtained from a direct 
numerical computation of 1^^ ° la^- 

5. Convergence of the series 

For T and / we now compute the radius of convergence of the expansion found. 
This can be done analytically. The series for the BirkhofF normal form and the 
symplectic invariant are then analysed numerically. 

To analyse the asymptotics of a„ define the ratio r„ = {2n + l)a„/((2n— l)a„_i 
and the recursion becomes r„ = (4 — n^'^)(K/2-|-r~_j^). The leading order iteration 
r„ = 2k + 4/r„_i has two fixed points at r = —2//) and r ~ 2p. The positive 
fixed point r — 2p is stable for positive k (or p > 1), while the negative fixed point 
r = —2/p is stable for negative k (or p < 1). The radius of convergence in h is 
given by l/|roo| which gives ^ min(/9, 1/p), and is thus controlled by the roots of 
ViP' closest to zero. For &„ a similar argument works after explicitly controlling the 
size of a„_i/6„_i and a„„2/^n-ii so both series with coefficients a„ and 6„ have 
the same radius of convergence. 



SEMI-GLOBAL SYMPLECTIC INVARIANTS OF THE EULER TOP 



16 




(a) Ratios tending to the (b) Radius of convergence as 

radius of convergence a function of k 



Figure 5.1. Radius of convergence of BirkhofT normal form and 
the symplectic invariant numerically. The results are the same as 
the data points overlap and curves coincide for both. In (A) we 
fixed K = 10 and varied N . In (B) we fixed N = 14, 15 in the ratio 
formula given in (A) and varied k instead, yielding the dark dashed 
(blue, N = 14) and light dotted (red, N = 15) curves respectively. 
The thick black curve is the graph of \ min(p, 1/ p). 



The Birkhofi^ normal form is given by the inverse of the series with coefficients 
dn/ {n + 1) which converges, we expect that this series converges as well. However, 
since inverting a series is a highly non-linear process explicit formulas cannot easily 
be obtained. Similarly, the symplectic invariant is based on the convergent series 
with coefficients 6„/(n + 1) — anjin + 1)^ composed with the Birkhoff normal form, 
so again we expect convergence. Numerically computing the ratio of successive 
coefficients as illustrated in Figure 5.1 indicates that the radius of convergence of 
both of these series is somewhat bigger than that of and the regular part of 1^. 
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